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ABSTRACT 

The transparent transport of real-time periodic traffic through 
a broadband packet network requires the recovery of source 
clock frequency at the destination. The inherently stochastic 
nature ofpacket transport in the network makes the recovery 
difficult. IT this paper, we illustrate the packet jitter 
phenomenon, describe several dock recovery algorithms and 
provide a systematic design procedure for the implementation 
of these algorithms. We apply them to a computer model of a 
single stage FCFS multiplexer to demonstrate their 
effectiveness and different design trade-offs. We also discuss 
the application of the algorithms to the audio and video 
services. 



1. INTRODUCTION 

The transparent transport of real-time periodic traffic across a 
broadband packet network requires synchronization of the 
network. A simple way to achieve this would be to distribute 
clock timing throughout the network via a common master 
clock. Alternatively, a less expensive approach would be to 
run the clocks at each node autonomously with a nominal 
frequency, and make the destination clock frequency follow 
the source clock frequencv as closely as possible. The source 
clock timing information can be explicitly earned in each 
packet header, as proposed for the narrowband integrated 
networks 131. for specific voice and video services on an end- 
to-end basis. The source clock timing information can also be 
derived from the arriving packet stream at the destination. 
This is achieved bv monitoring the destination buffer level, 
which is a linear function of the frequency difference between 
the source and destination clocks, and using it to adjust the 
destination clock frequencv without any buffer under- or 
overflow. These methods can also be combined together for 
some specific applications. 

Due to the statistical nature of packet transport in a 
broadband packet network, a periodic packet stream of a 
real-time source traffic tends to become quasi-periodic. The 
destination buffer level, therefore, cannot be directly used as 
a measure of the frequency difference between the source and 
destination clocks due to random fluctuations, which may be 
correlated as well as nonsialionary; its statistics may van- 
whenever there would be a change in the multiplexing and 
switching paitern due to addition or deletion of calls. 
The principal questions in the design of synchronization and 
clock recmerv schemes are: (i) how to separate the weak 
signal (the frequencv difference between the source and the 
destination clocks) from the large noise (the packet jitter) to 
obtain the estimate of the former; and (ii) how to control the 
destination clock frequencv using the frequency difference 
estimate. (The terms packet jitter and noise are used 
interchangeably in this paper.) The problem of 
svnchronizatioh and clock recover, in plesiochronous 
broadband packet networks has been discussed in |1.2| in 
rather general terms, with no specific design details. [4] 
addresses the problem in a quantitative manner. But their 
assumptions are quite restrictive, and the open loop control 
law suggested in [4] may lead to instability. In [5], a 



suboptimal and an optimal frequency difference estimation 
algorithms along with two control schemes for the adjustment 
of the destination clock frequency are described. These 
algorithms are based on realistic assumptions, and are 
different from those in [4]. 

In this paper, we focus our attention on the implementation 
of the alg orithms of [5], and present a systematic design 
procedure. Section 2 illustrates the phenomenon of packet 
Sitter via a simple example, where it is defined precisely in 
mathematical terms. Section 3 contains the estimation and 
control algorithms described in [5,6]. Implementation details 
are discussed in Section 4. Section 5 contains some 
simulation results by way of a computer model of a first 
come, first served (FCFS) multiplexer. Section 6 deals with 
service specific synchronization schemes, where the 
application of the above algorithms to audio and video 
services are discussed. Finally, Section 7 summarizes the 
results. 

2. PACKET JITTER PHENOMENON 
The phenomenon of packet jitter is illustrated in Fig. 1 for a 
single FCFS multiplexer with three packet streams of X,=1.0 
packet per second (mm). X;-2.0ppj and X,=4.0 ops (packets 
are assumed to be of fixed length) which are multiplexed to a 
a=8.0 ops line in the broadband network. The packets are 
removed svnchronouslv at the destination without the 
destination buffer over- or underflow at roughly the same rate 
as that of the source buffer, e.g., X,pjM , as shown in the 
figure. Fig. 2a and 2b are the plots of the destination buffer 
level as a function of time, which show that its value is zero 
at the indicated sampling instants when there is no jitter 
introduced into the third packet stream, and that its value 
deviates from zero at the sampling instants when the packet 
stream has to contend for emptv slots on the high speed line. 
(The positive slope in Fig. 2a is (ji-Xj) and the negative slope 
is X,.) 

Let / , be the source clock frequency in pps . (In the above 
example. /,-X,.) Let 6{k) be the destination buffer level 
function. The packets are removed from the destinauon 
buffer at the destination clock frequency of f : (k)pps. Let 
^/ (*)=/,-/,(*) be the difference between the source and 
destinauon clock frequencies. The rate of change of the 
destination buffer level indicates the measure of Af [k). 
When the packets in the input line arrive with sti>chastic 
delavs. the input packet stream represents a noisy value of 
the source clock frequencv, /,. and o(* ) fluctuates in a 
random manner depending on A/ (A) and the packet jitter 
characteristics. Let the noisy frequency be f,+n(k). where 
n (k ) is the random fluctuation in the source clock frequency 
due to the the stochastic inter packet spacing. Then the 
frequencv difference between the incoming and the outgoing 
packet streams to and from the destination buffer is 
V (fcl=A/(* )+"(*)• The dynamics of the state 

equation for o(k ) can now be derived from the conservation 
law. A convenient choice of the sampling interval is 
7"(* )=l//^(k). the time period of packet removal from the 
buffer. 

o(k + 1 )=o(* )+ T(t )(/ 1-/ :(*)+"(* ) (: - 1 > 
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i*(tH(t+iM*)=''(*)A/(tM(*) (") 

where d (k )=n (k )T(k). d(k) is the equivalent noise term in 
number of packets. The nonlineanty arising from the choice 
of T(k ) is ignored. (2.1) or (2.2) is the state equation of the 
buffer. 

Definition of packet jiaer: We define the random variable 
d(k) as the packet jitter. When ^/(*)-0, the random 
change in 4(-) from one time instant to another represents the 
packet jitter. 

The most difficult part in the design and analysis of the dock 
recovery schemes is the modeling of the noise process, d(k). 
Its expected value, of course, is zero, i.e. Edlky^O (£• is the 
expectation operator.), by the fact that if the frequency 
difference is zero, then the average accumulation of packets 
in the destination buffer must be zero, assuming no packet 
loss through the network. The other statistics, in general, are 
not known, however. Some preliminary analysis indicates that 
d(k)'a expected to be a correlated sequence. The correlation 
is dependent on several factors - loading of the network; 
relative magnitudes of \ t 's in comparison to n; service 
disciplines at the multiplexers and switches, e.g., FCFS or 
train scheduling [11; type of connection, e.g., virtual circuit or 
datagram. These factors will influence the correlation of the 
packet jitter in a complex manner. However, in general the 
correlation in the jitter sequence results from a repeating 
pattern of jitter every few packets in the stream. If the 
packet stream goes through manv stages of multiplexing and 
switching, the period of the pattern is large. This would mean 
less correlation. The plot of Fig. 2b shows this repeating 
pattern of packet jitter every 4 packets for the example 
considered above. 

For simplicity, in this paper we do not explicitly consider the 
correlation of the process d (k) in design. The jitter process, 
d(k), is assumed to be stationary with Ed (k )=0 and an 
unknown constant variance, denoted by v . 



3. CLOCK RECOVERY ALGORITHMS 

We assume that the network is plesiochronous. i.e., the 
clocks at the source and the destination derive their timings 
from independent local clocks, with the same nominal 
frequency but different ppm values. The frequency difference 
between the two end clocks. A/(t). is quite small - typically 
a few tenth of a ppm . But the values of d(k ) may vary from 0 
packet to a few packets from one instant to another. Thus, 
the changes in o(k) are mainlv due to the random jitter d(k). 
However, since £d(fc)=0. we expect that the contribution of 
the ±f(k) term to the buffer level function in j()T(k) time 
periods mav become quite significant, where /(•) denotes the 
length of the interval, which is large and time-varying. Thus, 
it mav be possible to estimate the frequency difference from 
the data set o\k). k>0 if we wait long enough, i.e., for long 
enough / (-). 

We propose the following two time-scale model - a fast time 
sequence, k >0. when we collect observation data o{k ), and a 
slow time sequence, m >0, on which we compute, the 
estimates of the frequency difference, denoted by A/(m). 
from the set o(*). *>0". and adjust /;(m) using this 
estimated frequency difference. The slow time scale interval 
|m . m + l) is made timc-varving. i.e.. (KKj(m) for the slow 
lime scale interval. We propose to use a lime-varying lime 
period governed by 

;(m + l)=-^y;(m) < 31 > 

where i(m) is the length of the interval [m , m + l), with /(0) 
as its initial value, (l-o) is a design parameter which 
controls the rate of decay of Af (m) as described later. The 
size of the interval increases at the same rate at which ±j (m) 
decreases. The idea behind using this scheme is that in the 
beginning when the frequency difference is large, we want 
more frequent control actions, perhaps at the cost of poor 
estimate. As time progresses, the frequency difference 
diminishes exponentially at (l-o) rate, but the noise acts 



persistently, resulting in a decreasing signal to noise ratio with 
time. Therefore, we must have a larger interval for a better 
estimate of the decreasing A/ (m ) with m . 

We process the observation set *(*), <X*<y'(m) in some 
manner to filter out the effects of the rapidly changing (from 
packet to packet! noise and obtain the frequency difference 
estimate a/ (m ) for each interval [m, m+l), which is used in 
the control algorithms below to compute / 2 (m+l). Note that ; 
±f(m) is constant in each interval [m,m+l). Let A/(m); 
denote the frequency difference estimation error, i.e., : 
i/(ra)» i/(mj-i/jm). The estimation algorithms 
described later are such that the estimate A/(m) is unbiased, 
i.e., £A/(/n)-£A/(m). Further, if d(k) is a white noise 
sequence, then the resulting A/(m) is also a white noise 
sequence. 

3.1 Frequency Adjustment Algorithm 

Before we describe the algorithms to estimate the frequency 
difference, let us discuss the control algorithms. Since 
Sf (m ) is not known, its estimate A/(m ) is used as an input 
to the control algorithms. 

Control Law I: A simple control law is to add a weighted 
value of A/(m ) to J 2 {m ) to obtain / 2 (/n +1), which results in 
a first order dynamic 'system, i.e., 

/,(m+l) = / : (m) + oA/(m) (3.2) 

for some 0 < a < 1. With the above definition of A/(m), 
the closed loop system is 

/.(m + l) = (l-tO/jfmKo^+A/fm)) 

which is stable for the choice of a. Since £/i=/i and 
£ A/(m M>, the mean of f z (m ) satisfies 

£/ ; (m + l) = (l-a)£/ : (m) + o/, (3-3) 

whose asymptotic value can be shown to be £/ 2 (oc)=/,. 

The dynamics of the variance »(m) = E{JA.m)-EJ z (m)f 
can be' derived from (3.2) assuming that ^/(m) is a white 
noise sequence, and is given by 

r(m + l) = (l-o) : x(m) + o ; <t("») (3.4) 

where <r(m )=E±p(m) is the variance of the frequency 
difference estimation error. Since 0 < (l-o) < 1, the 
asymptotic variance (as m — »oc) of f z (m ) can be shown to be 
*(x) = <r(x). c{x) is related to v through the 

estimatiorTalgoriihms described later. The residual variance 
j(x) is smaller than a(x) for all 0 < a < 1. There is a 
trade-off between how fast the dynamic response should be 
and him much variance reduction' we can obtain. A larger n 
means faster response ((l-o) is closer to zero), but higher 
residual variance, and vice versa. Also, a larger a implies a 
smaller destination buffer. 

Control Law II: The control algorithm, in which the value of 
/-■(m + l) is obtained bv adding the weighted values of / ; (m), 
/^(m-1), A/ (m) and A/ (m— 1), provides more freedom in 
controlling the rate of convergence and the asymptotic 
variance. The weightings on / j(m) and ( 2 (m-l) should be 
such that the open loop transfer function should contain 
(l_r-i) as a factor in the denominator, so that there is no 
steadv state error, (z is a forward shift operator, i.e.. 
w (f )=z(/+l).l The new control algorithm results in a 
second order dynamic system 

/ 2 (m+l) = (l+i)f z (m)-6f ,(m-l)+a±f (m)+J±f (m-l) 

where 4, a and 3 are free parameters, which can be 
determined on the basis of the choice of the two closed loop 
poles, while minimizing »(oe). See (5,6] and the references 
therein for detail. 

3.2 Frequency Estimation Algorithms 
We describe below two increasingly complex algorithms for 



computing <l/(/n). Both algonthms provide «rimates °f 
l/(E)by processing the finite dau set A*(*)i,0$*^">2 
me interval [m,m+l). For convenience, theTjuffer state 
equation (2.2) is rewritten here as 

Simple Ttoe-ATeragta* Estimation Algorithm: All ad hoc 
scheme to estimating A? (m^ may be just to take the sample 
mean of all the available observations ^(■)f 1 (m). This is 
motivated by the idea that the zero mean noise samples will 
cancel out by averaging it. There are a total of y(m) samples 
in the interval [m , m+1), and thus, 

IK" 1 ) kmO 

./M r(ffl+ l)-*(m)) P-6) 
/(«) 

where <r(m+l) is the value of «(•) just before the control 
action is taken at the instant m+1. 

Optimal Estimation Algorithm: The optimal estimation 
algorithm makes use of t(0) and v. and is based on the first 
control law. The Appendix contains a mechanism by which 
the variance of the noise is estimated in real time using the 
same set of data 4d(). The following set of equations for 
the optimal estimation algorithm are developed in 15.6] from 
the basic Kalman filter equations and the first control law 

*/(m ) P , < , "1 ) , . («-(w + l)-«(m ))/;(mm 

,y ' l+i(m)p(m) 

+ -Elfim) 

,(m + l) = O-^l^ 1 ,^^, 



(37) 
(3.8) 



with p(0)= 



/j(0)v 



known. d~(m+l) is the value of the buffer 



1 function at the instant just before m+1. £A/(m), /(m) 
/ 2 (m) are all known quantities in [m, m + 1). For the 



level I 

purpose of computation, we approximate f,(m )-f i(0). 
Since v is not known, we propose to use the variance 
estimation algorithm, (A.l) or (A.2), to estimate it in the 
interval (0,1) (assuming stationariry), which wul be used 
during the other intervals as well. (3.9) is obtained from the 
quadratic Riccati equation of the Kalman filter equations (see 
[5,61 for detail.). 

In both cases, the algorithms result in finite estimates of 
frequency difference tor each interval, and <x(m)— 0 as 
i(m)—<x> with m— oo if d(k) is assumed to be a white noise 
sequence. The first algorithm is applicable to all systems 
regardless of whether d(k) is correlated or not. The second 
algorithm may give poor performance if the sequence d(k ) is 
hShlv correlated. As said earlier, modeling of correlation of 
5ft)' is difficult. We are continuing our efforts to understand 
the various contention mechanisms which cause correlation in 
packet jitter. Even then, it does not appear straight forward 
to develop optimal estimation algorithms for the correlated 
<f(t)case. 

Fig. 3 shows the signal flow diagram for the estimation and 
control algorithms along with the two time scale model. 

4. IMPLEMENTATION OF THE ALGORITHMS 

Let /,=/±/c,, where / is the nominaj frequency and ^, is 
its standard deviation in ppm , i.e., / <ft<f + "^» ' 

Likewise, let f < /:(<>) < Z+^TT • TOs su S8 es,s ,nal 

f 

Efi m £/ : (0) - / and the variance of / , is and that of 



/ 2 (0) is ^jjj-. Therefore, £4/(0)-0 and i<0) ^5— 

should be used in the optimal estimation algorithm, once it is 
ascertained that d(k) is close to being a white noise 
sequence. 

Now. the next question is how to choose ;"(0)? i'(0) is 
determined from the prior knowledge of 4/(0) (or w(Q)), the 
maximum value of the noise magnitude and the allowable 
buffer size, so that even if the system is running for ; (0l units 
of time in the presence of noise, there is no buffer under- or 
overflow, and we get a reasonable estimate of 4/(0) at the 
end of the first interval. In practice, after a few iterations 
A/(m) will become verv small. At this time we can stop 
changing / (m ) and use a constant value. We propose that for 
some m-m,, whenever either |.i/(/n,) | «„ riWi)* 4 !- or 

* then set /(m )-/("> + 0 for m>m x for the 

v(mi) . . 

steady sute operation. Where <„« 2 and «, are design 

parameters. 

In the above algorithms, the time intervals [m, ffl+1). m>0 
are determined a priori, and the destination buffer levels are 
read at time instants, m >0. The values of these may contain 
a fraction of packets. Alternatively, it is also possible to fix 
the destination buffer level function at discrete packet levels, 
e g half, three quarter etc.. and to measure the time 
intervals for o(k) to attain these predetermined labeled 
values. There are two limitations of this technique: (i) at*) 
mav occasionally cross the labeled points due to a large but 
unexpected jitter. If for some reason these crossings occur 
too frequently, the estimation of the frequency difference may 
be quite poor due to small number of observations, much 
will result in a larger residual jitter, (u) By observing *(* ) «» 
a coarse level, we loose information about the statistics of the 
high frequency jitter component buried in it. Thus, the 
vananceof the jitter cannot be estimated on-line from this 
observation set, and the Kalman algorithm cannot be used. 

5. SIMULATION RESULTS 
We used a computer model of a FCFS multiplexer to study 
the estimation and control algorithms of Section 2. lac 
inputs to the multiplexer consist of a combination of signals, 
including DSO andDS3. The output rate of the multiplexer is 
150 Mb ps. We implemented both estimation algonthms ((3.6) 
as well as (3.7) and (3.8)) in the first control algorithm (3.2) 
with the parameter Aarytag between 0 and 1 to ^^.^j 
trade off between the speed of convergence and the residual 
variance. The algorithms are implemented to synchronize 
clocks at DSO and DS3 rates - the two extremes of the DS- 
hierarchy. The various values used in the simulations are: 
For DSO signal - / ,-64.00 kbps» /,(0)-63.99 kbps and 
55 (OpMwn -, «nd tot DS3 signal - / ,-44736.00 kbps. 
/ j(0)-44735.11 kbps, and 4/ (0)-20 ppm . 
The choice of DSO and DS3 signals was dictated by the 
expectation that they generate jitter with different degrees of 
correlatedness. We first compute the autocorrelation function 
of the jitter generated by these two signals. Figs. 4 and 6, 
respectively, show the plots of normalized a"^""* 13 "™ 
function (normalized by variance) for the DSO and DS3 
sumals. We observe that the magnitude of the jitter 
generated by the DSO signal is small (standard deviation - 
0 006 packets) with a decaying correlation. (See Fig. ■».) me 



signal to noise ratio is 



=0.026. On the other hand. 



the DS3 signal provides^ C |arge jitter (standard deviation - 
0 26 packets) with a large and persistent correlation and a 
shorter period of jitter pattern. (See Fig. 6.) The signal to 
noise ratio in this case is °° 00008. which is 

much smaller than before. 

As expected, therefore, the estimation algorithm (3.7) with 
the control law (3.2) results in a superior performance (see 
Fig. 5a) in terms of reduced jitter for the DSO signal 1 as 
compared to the ad hoc scheme (3.6) with (3.2) (see rig. Sby 
The value of o-0.5. Solid lines denote the estimated 




frequencies and dashed lines represent the buffer levels. 

Figs. 7a, 7b and 7c show the plot of the first estimation (3.6) 
and control algorithm (3.2) for OS3 signal with different 
values of a ; a - 0.25, 0.50, and 0.75. In these plots, we 
observe the effect of the parameter o on the convergence 
speed and the rapid fluctuation in -l/(m). With a larger a, 
tne estimated frequencies converge to the source frequencies 
faster but show greater fluctuations. However, because the 
estimation and control algorithms correct the frequency 
difference quickly when a is large, the system operates at 
smaller buffer levels despite the variations in the estimated 
frequency. (See Fig. 8j Since the packet jitter is highly 
correlated for the DS3 signal, the optimal estimation 
algorithm results in very poor estimates, and consequently 
poor convergence. 

6. SERVICE SYNCHRONIZATION 



The clock recovery schemes presented earlier do not utilize 
the characteristics of individual services such as voice or 
video. Moreover, they are sensitive to the packet loss and 
large stochastic network delays. Large network delays mean a 
larger fixed delay for the original signal. A loss of packet 
means that the incoming packet stream is not a true measure 
of the source clock frequency. Thus, it appears that the 
potential application of the above clock recovery schemes will 
be limited to the packet transport of the existing digital 
circuit-switched hierarchical trunk signals. For an end-to-end 
packet-switched network, one can employ more flexible 
service-dependent clock recovery schemes, utilizing the 
properties of the individual services. 

6.1 Voice Signal 

The end-to-end voice synchronization scheme presented here 
is different than the ones in the literature [3]. For the 
synchronization of a voice signal, it is not necessary to use 
tne clock recovery schemes of the previous section to adjust 
the destination clock frequency. The destination buffer under- 
and overflow due to the frequency difference between the 
source and the destination clocks can be controlled by 
occasionally deleting or repeating packets in the voice packet 
stream at the receiving end after delaying the packet stream 
by the maximum jitter amount. For example, with a clock 
frequency difference of 10 ppm , the average service 
interruption interval (average period of packet deletion or 
addition) would be more than 26 minutes with a 1000 bit 
packet length for 64 kbps voice call in a 150 Mbps packet 
network. Moreover, in practice, this clock difference can be 
corrected bv deleting or adding packets during the silence 
interval without affecting a normal telephone conversation. 

6.2 Video Synchronization 

Depending on the video signal rate, type of coding and the 
desired quality of the video service, one may use any of the 
following methods for synchronization " of end-to-end 
terminals. 

6.2.1 Without Clock Recovery Schemes: 

Similar to the previous case, independent local clocks with 
the same nominal frequency can be used to achieve 
synchronization between a video encoder at the source and a 
decoder at the destination. The incoming packets can be 
delaved bv a given amount of time (more than the maximum 
tolerable jitter, possibly one video frame delay) and played 
out according to the local decoding clock. Terminal 
synchronization can be achieved by repeating or deleting an 
entire video frame at the destination buffer, depending on the 
buffer level. This is like the phenomenon of slip in the digital 
circuit switch interface. If the intra-frame coding is used, 
then packets corresponding to one whole frame can be 
deleted or repeated prior to decoding. This can be 
accomplished by inserting special frame sync pattern in the 
frame boundary at the transmitter to distinguish one frame 
from the other.' If the clock difference is I ppm . then frame 
repetition or deletion will happen once in about even- 9 hour 
interval, which is acceptable for most of the video services. 



6.2.2. With Clock Recovery Schemes: Fixed Bit Rate Video: 

Straightforward PCM or fixed length DPCM type of simple 
coding schemes generate a constant bit rate, and 
packetization of this coded bit stream results in a periodic 
packet stream at the source. The clock recovery scheme 
described in the previous section can be effectively used to 
synchronize clocks at the two ends, assuming a reasonable 
jitter and no packet loss. 

Variable Bit Rate Video: Variable bit rate generation is 
inevitable for almost all video ceding algorithms that utilize 
local image statistics, except for straightforward PCM or 
some low compression coding. This coding scheme seems to 
be a natural choice for a packet-switched network. In this 
case, however, the clock recovery schemes of the previous 
section cannot be used directly with the variable rate packet 
stream, since it is difficult to distinguish the variation in the 
periodicity of the packet stream due to the statistical coding 
and due to the stochastic network delays. Therefore, the 
source has to provide an extra timing information to the 
destination, so that the clock recovery circuit can lock on to 
that timing information. The timing information can include 
sequence number (video line or frame number) and a special 
pattern to easily detect the sequence. This tuning information 
can be inserted into the variable bit stream prior to 
packetization. If the packet length is short, it is also possible 
to insert separate timing packets into the packet stream, and 
the source clock frequency can be recovered from this 
jittered timing packets based on the schemes described in the 
previous section. 

Subband Coding: In the subband coding scheme an input 
video signal is decomposed into several different frequency 
bands, and different coding schemes are applied in each 
frequency band to take advantage of the local image statistics. 
This generates multiple parallel bit streams. The idea behind 
the clock recovery for such video signals is to apply a fixed 
bit rate coding scheme in one band, possibly the lowest 
frequency band, and recover the source clock frequency from 
the packet stream coming from that band at the destination. 

7. CONCLUDING REMARKS 

The difficulty in clock recovery for real-time periodic traffic 
in a broadband packet network stems from the random 
nature of packet transport through the network. We have 
presented several schemes by which the stochastic effects on 
the periodic packet stream can be mitigated, and the source 
clock frequency can be recovered asymptotically. We have 
discussed the application of these schemes to the audio and 
video services, as well as to the digital circuit-switched traffic, 
e.g., DS1, DS3. A systematic implementation procedure has 
been presented, and' the effectiveness of the schemes have 
been evaluated via simulations bv way of a computer model 
of a single stage FCFS multiplexer for different source traffic 
rates. 
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APPENDIX 

On-line Computation of Noise Variance Estimate, £(-)-' v(k) 
can be estimated from the available information, 
M(k +1), A*(k ) Ad(0) as follows. , 

Approximate Analysis: Since A/ (k) is small, we can ignore j 
its effect on Ad(fc) on the fast time scale and write 
A*(t) = d(i) for all *>0. Hence, v - Ed\k) = EA* 2 (Jt) ; 
Now, the variance, v, can be estimated as follows (v(k) is • 
the estimate of v at time k .) 



- *tr ' ( * } + rh^** (A1) 

for all oo>*>0, with ■7(0>»A4 2 (0), or some initial guess of 
0(0), if it is available. 

Exact Analysis: Successive differencing of the observation 
Ac*( ) yields AA<*(t) > A4(*+l)-A«(*)-<f(ft+l)-<t (*). 
The expected value of AA<4(*) is zero, and its variance is 
given bv £A^*)-£(rf(*44)-rf(*))*- Since the noise is 
uncorrelated and stationary, this yields 
EdHk+ll+EdHk)* El&4>\k)=2y. Analogous to (A.12), 
the variance estimate, v, can be computed via the following 
recursion 

Hk+l) ° kTi Hk) + iir 2^^ k) (A2) 

with i?(0)= ^ as the initial condition. 
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